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DETERMINING MODES FOR THE 
3D NAVIER-STOKES EQUATIONS 

ALEXEY CHESKIDOV, MIMI DAI, AND LANDON KAVLIE 


Abstract. We introduce a determining wavenumber for the forced 3D Navier- 
Stokes equations (NSE) defined for each individual solution. Even though this 
wavenumber blows up if the solution blows up, its time average is uniformly 
bounded for all solutions on the global attractor. The bound is compared to 
Kolmogorov’s dissipation wavenumber and the Grashof constant. 
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1. Introduction 


The Navier-Stokes equations (NSE) on a torus, are given by 


, , f u t + (u ■ V)m - uAu + Vp = f 

1 j \ V • it = 0, 

where u is the velocity, p is the pressure, and / is the external force. We assume 
that / has zero mean, and consider zero mean solutions. 

The dissipative nature of these equations is reflected in the existence of an ab¬ 
sorbing ball in L 2 . Moreover, in the two-dimensional case, there exists a compact 
global attractor which uniformly attracts all bounded subsets of L 2 . This attracting 
set is, in fact, finite dimensional, as was proven by Foias and Temam in |14j (see 
also Cl)- The first result for the finite dimensionality of a two-dimensional fluid 
appeared in the work of Foias and Prodi m, where they showed that high modes 
of a solution are controlled by low modes asymptotically as time goes to infinity. 
The number of these low modes, called determining modes, was estimated by Foias, 
Manley, Temam and Treve [T2] and later improved by Jones and Titi [20] . See also 
maim nz] and references therein for related results. 

In three dimensions the situation is drastically different as the equations have 
thus far eluded a proof for the existence of classical solutions. Even so, the existence 
of a global attractor for weak solutions is known in a weak sense mm- This weak 
global attractor consists of points on complete bounded trajectories and attracts all 
bounded subsets of L 2 in the weak topology. However, it is not known whether the 
solutions on the attractor are regular, unless the attractor consists of a single fixed 
point. Neither is it known whether the attractor is compact or finite-dimensional. 
Similarly, the existence of a finite number of determining modes is not known in the 
three-dimensional case. Nevertheless, Constantin, Foias, Manley, and Temam [6] 
showed the existence of determining modes assuming that the H 1 norm of solutions 
is uniformly bounded. The question whether the global attractor of the 3D NSE 
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is bounded in H 1 is open and may very well require a resolution of the regularity 
problem. However, even assuming regularity, this would not immediately guarantee 
that the H 1 bound would depend only on the size of the force (a Grashof constant), 
and not on the shape of the force. 

With no hope of getting a finite number of determining modes for the 3D NSE, 
one might ask whether this can be done in some average sense. Indeed, the Kol¬ 
mogorov 41 phenomenological theory of turbulence [21] predicts that the number 
of degrees of freedom should be of order where Kd is Kolmogorov’s dissipation 
wavenumber. This number is often used as the resolution needed for direct numer¬ 
ical simulations, so one might ask an alternative question: What is the number of 
determining modes for a time discretization of the 3D NSE and how does it depend 
on the force and time step? 

In this paper, without making any assumptions regarding regularity properties 
of solutions or bounds on the global attractor, we prove the existence of a time- 
dependent determining wavenumber A u (t) defined for each individual solution u. 
We show that any weak solutions on the global attractor u and v that coincide below 
max{yl„, A v } have to be identical. The wavenumber A u (t) blows up if and only if 
the solution u(t) blows up. Nevertheless, the time average of this wavenumber 
is uniformly bounded on the global attractor, which we estimate in terms of the 
Kolmogorov dissipation number and Grashof constant. 

To begin, let u be a weak solution of the 3D Navier-Stokes equations. For 
r £ (2,3) we define a local determining wavenumber 

A u , r (t) := min{A, : A p r ||u p ||ia < cvi', \/p > q and X~ ||u< q ||L~ < c r v , q € N}, 

where c r is an adimensional constant that depends only on r. Here \ q = ^ with 
L being the size of the torus, and u q = A q u is the Littlewood-Paley projection of 
u (see Section [2]). Thanks to Bernstein’s inequality, we have 

A u ,r > A^ ls := min{A ? : A“ 1 ||u p ||i,«. < cqu, \/p > q, q £ N}, 

which is a local dissipation wavenumber introduced by Cheskidov and Shvydkoy in 
[8] . It defines a dissipation range where a local Reynolds number corresponding to 
high frequencies is small, i.e., 

K := gKiL °° < co, VA 9 > , 

H V 

where l q = A” 1 . The dominance of the dissipation term above is reflected 
in improved Beale-Kato-Majda and Prodi-Serrin criteria where u is replaced with 
its projection on modes below yl)( ls (see 0). The determining wavenumber A u>r 
imposes tougher condition on high modes, as well as requires a control on low 
modes via the low frequency Reynolds number 


K ■■= 


q \\U< q \\ L ° 


< -C C-f 


X q — A ur . 


It is also worth mentioning that a similar determining wavenumber is used in [1] 
to prove the existence of a finite number of determining modes for the surface 
quasi-geostrophic equation equation in critical and subcritical cases. Even though 
the determining wavenumber enjoys uniform bounds in those cases, it still proved 
useful to start with a time dependent wavenumber defined based on the structure 
of the equation only, and then study its dependence on the force using available 
bounds for the global attractor. 
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We prove the following. 

Theorem 1.1. Let u(t ) and v(t) be weak solutions of the 3D Navier-Stokes equa¬ 
tions. Let Aft ) := max{yl Ujr (t), A Vtr (t)} for some r £ (2,3). Let Q(t) be such that 
A{t) = A Q (t) . If 

(1-2) u(t)<Q(t) = v(t)< Q{t) , Vf > 0, 

then 

lim ||u(t) - v{t) \\ L 2 = 0. 
t—yoo 

We prove this theorem in Section [3] A minor modification of this result gives us 
the following theorem: 

Theorem 1.2. If u(t) and v(t) are two Leray-Hopf solutions on the weak global 
attractor A such that 

(1-3) u(t)<Q(t) = v(t)< Q(t) , Vi< 0, 

where Q is given in Theorem \ 1. 11 then 

u(t) = v(t), Vf < 0. 


This establishes the existence of determining modes for the 3D Navier-Stokes 
equations. It is worthwhile to note that the determining wavenumber A U}V depends 
on time and may not be bounded. Actually, it is bounded if and only if u is 
regular. However, the average determining wavenumber (A) = A J) t+T A u , r (r) dr 
always enjoys a uniform bound. Indeed, we establish the following pointwise bound: 


(1.4) 


Au.rft) 


< 


l|V«(t)|| 


L 2 


Note that this automatically provides a finite number of determining modes and 
recovers the results by Constantin, Foias, Manley, and Temam [5] in the case where 
l|v«(i )||£ 2 is bounded on the global attractor, which is known for small forces 
(laminar regimes). On the other hand, (11.41) holds in general for arbitrary forces 
and implies that (A) is uniformly bounded for all Leray-Hopf solutions on the global 
attractor, i.e., complete bounded trajectories. However, the bound m is sharp 
only in the case of extreme intermittency, where on average there is only one eddy 
at each dyadic scale. To make a connection with Kolmogorov’s turbulence theory 
ED DU we have to define an intermittency dimension and analyze (A) in various 
intermittency regimes. 

In Section [4] and in Section [SJ we further examine (A), comparing it to Kol¬ 
mogorov’s dissipation wavenumber as well as the Grashof constant, defined as 


Kd := 


(£)’ 


G := 


I H- 


aa 


1/2 


e:=\ d 0 v{\\Vu\\l 2 ), 


where d £ [0, 3] is the intermittency dimension.This parameter is defined in Sec¬ 
tion U in terms of the level of saturation of Bernstein’s inequality. The case d = 3, 
where there is no intermittency and eddies occupy the whole region, corresponds 
to Kolmogorov’s regime. In this case the bounds read 

4+5 


(^>< 


5 


2+<5 \ -1-5 


{A) A 0 


G 


(tva§ 


+ G" 


d = 3, 
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where 8 can be arbitrary small when r is chosen close to 3. On the other hand, in 
the case of extreme intermittency, the bounds become 

(A) iS K d, (A) 1 $ „ 9 ,—F AoG 2 , d = 0 . 

TV Ao 

Finally, in Section [ 6 ] we show that if one of the weak solutions in Theorems ll.il 
fi~2l satisfies the energy equality, then A can be defined in terms of that solution 
only. The energy equality holds for regular solutions, such as steady states, or 
solutions belonging to Onsager’s space L 3 (0, oo; BlJ^) (see [3] where the proof can 
be applied to the viscous case, where the space B-J^ can be replaced with B 3 ^). 


2. Preliminaries 


2.1. Notation. We denote by A < B an estimate of the form A < CB with some 
absolute constant C, by A ~ B an estimate of the form C\B < A < C 2 B with 
some absolute constants Ci, C 2 , and by A < r B an estimate of the form A < C r B 
with some adimentional constant C r that depends only on the parameter r. We 
write II ■ lip = II • IUp, and (■, •) stands for the L 2 -inner product. 


2.2. Littlewood-Paley decomposition. The techniques presented in this paper 
rely strongly on the Littlewood-Paley decomposition. Thus we here recall the 
Littlewood-Paley decomposition theory briefly. For a more detailed description 
on this theory we refer the readers to the books by Bahouri, Chemin and Danchin 
PQ and Grafakos fl9l . 

We denote X q = for integers q. A nonnegative radial function \ G c5?°(R 3 ) is 
chosen such that 


(2.5) 



for iei < | 

for |f| > 1 . 


Let 

<^(C) == x(C/2) - x(0 


and 

for q > 0 , 

X(C) for q = -1, 



so that the sequence of ip q forms a dyadic partition of unity. Given a tempered dis¬ 
tribution vector field u on T 3 = [0,L ] 3 and q > —1, an integer, the qth Littlewood- 
Paley projection of u is given by 


u q (x) := A q u{x) := ^ u{k)4> q (k)e l2 Z k ' x , 
fee z 3 


where u(k) is the fcth Fourier coefficient of u. Note that it_i = u(0). Then 


OO 

u = J2 u i 

q =-1 


in the distributional sense. Note that 


INI H* 



1/2 
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for each u € H s and s£l. To simplify the notation, we denote 

Q 

u<q := u q- 

g= —1 

2.3. Bernstein’s inequality and Bony’s paraproduct. Here we recall useful 
properties for the dyadic blocks of the Littlcwood-Paley decomposition. The first 
one is the following inequality: 

Lemma 2.1. (Bernstein’s inequality) Let n be the spacial dimension and r > s > 1. 
Then for all tempered distributions u, 

( 2 - 6 ) IKIk<A? ( —Vjs- 

Secondly, we will use Bony’s paraproduct formula. First, note that 

U • V V — ^ ^ U<(p —2 * ^ Vp ^ ^ Up • V V<fp —2 + E E Up ■ Vty . 

p p v b-p'|<i 

Due to (12.51) we have = 0 when |£| < 3/4 or |£| > 2, and hence 

ifq9<q—2)>q+2 = 0, (/<?3<g—2)<g—3 = 0, <7qr+i) ><j+3 = 0, 

for tempered distributions / and g. Therefore, 

A q (u-X/v) = ^ A ? (m < p _2 ■ Vn p ) + A g (rip ■ Vu<p_ 2 ) 

g—l<p<g+2 g—l<p<q+2 

+ 'y ) ^ ( Ag (Up ■ Vfp' ). 

p>q-2 |p-p'|<l 
p>q-2 

It is usually sufficient to use a weaker form of this formula: 

A, (u-Vv)= ^2 Ag(w<p_ 2 • Vup) + y^ Ag(rtp • Vu<p_ 2 ) 

I«-pI<2 I«-pI<2 

+ ^ ' Ag({(p.VUp), 

p>q-2 

where ftp := it p _i + u p + it p +i- 


3. Proof of Theorem o 
Let w := u — v, which satisfies the equation 
(3.7) uit + u ■ Vu> + w ■ Vv = — V7r + vAw 

in the distribution sense. By our assumption, ru<Q(t) ( t ) = 0. 

Multiplying (13.71) by Agi/;, integrating, and adding up for all q > —1 yields 

ylM 0 ll 2 - ^IKto)||2 + u II Vrc || 2 dr 

(3 8) < / > / A q (w • S7v)w q dx dr — / / / A q (u ■ Vw)u> q dx dr, 

Jto J t 3 Jt 0 ,77^ J T 3 


?>-i 


9>-l' 


=: / /i dr + 
Jto 


h dr. 


to 
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Using Bony’s paraproduct mentioned in Subsection 12.31 I± is decomposed as 


h = - Y, Y, I A q (w<p -2 ■ Vv p )w q d,X 

9>-l |g-p|<2' 


/ T 3 


- "y 'y I A q (w p ■ S7v<p_2)w q dx 

q>-l \q—p\<2 ' "" 


/T 3 


-EEL A q (ui p ■ \7v p )w q dx 

q>-l p>q-2 ' 


/ T 3 


—:Jn + I \2 + 1 13 . 


These terms are estimated as follows. Recall that w<Q(t) = 0. Then using Holder’s 
inequality we obtain 

ms E E / \Aq(w < p -2 ■ Wv p )Wq\dx 

p>Q+3\q-p\<2 JT3 

< E E ll w <P-2||^Ap||i;p|| r ||u; 9 ||2. 

p>Q+3 \q—p\<2 

Now using definition of A VtT , (12.61) . and Jensen’s inequality, we get 

\hi\< c rv y y y~~\\wqh y iiuyii^ 

p>Q+3 \q— p|<2 p'<p—2 

<c r v y A ? ||w 9 || 2 y Ap/||wp/|| 2 AjZJ, 

q>Q+l p‘<q 

< c r is\\\7w\\l, 

where we used the fact that r < 3. 

Next, we have 

i /i2 i -EE / |A q{Wp ■ Vv<p- 2 )Wq\dx 

p>Q\q-p\<2 JT3 

^E E / |A g (w p • Vv<Q)w q \ dx + E E E /J A q (Wp ■ S7Vp')w q \ dx 

p>Q\q—p\<2 T3 p’>Q p>p'+2 \q—p\<2 T3 

< iiv^qHoo y y ikihkii.+ EiivvIIoo y y ikimkiu 

q>Q \q—p\<2 P’>Q p>p'+2 \q—p\<2 

< iiv«< Q iioo y \\w q \\i + y Ap^iMir y ikh 2 

q>Q p‘>Q q>p’ 

Using definition of A Vjr and Jensen’s inequality, we obtain 

i j i 2 i < c r v y aqa q 11 Wq hi+ c r v y a i> y \\w q \\i 

q>Q p’>Q q>p ' 

<c r ^A 2 |Kii 2 + C r^A^iKiii y \-y 

q>Q q>Q Q<p'<q 

< CrI/||VlI7||l. 
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We will now estimate / 13 . From w<q = 0 and w<q- 1 = 0, it follows 

M ^ E E / \A q (w p ® f P )V'U) g | dx 

q>Q p>q—2 T3 

= ^2 ^2 / lAgC^P <8>u p )V«; 9 |da; 

p>Q Q<9<P+2 ^ t3 

= E E J A g (w p < 8 > u p )Vw g | dx + E /JA 9 (w}q (8 WQ)Vw; g | dx 

«/T 3 «/T 3 


p>Q Q<q<p +2 u Q<g<Q +2 ' 

^ E H^pIMMIr E A 9 ||w g ||_^i_ + Xq ^2 11 112 11 XQ 11 00 11 Rig 112 


p>Q Q<q<p+2 

We use definition of A lsr to conclude 

, i-5 11 ~ 11 , i+5 


Q<q<Q +2 


I^13| < CrV A p r ||wp|| 2 X! A 9 +r ||w 9 || 2 + C r j/A^ ^ II^Qlblkglb 

p>Q Q<«?<P+2 Q<g<Q+2 


< Cr ^ApK|| 2 ^ A g ||w g || 2 A g r _ p + c r v\ 2 q 


\w. 


q 112 


p>Q 


Q<q<p +2 


Q<9<Q+2 


< 


Crl/IIV^Ha. 

Therefore, we have 

(3.9) |/i|<Cri/||V«;||i 

To estimate / 2 , we start with the decomposition 


h = ~ E E / A q (u<p-2-Vw P )w q dx 

q >-1 k-p|<2‘ " Q 


/T 3 


-EEL A g (iip • Vui<p_ 2 )w g da; 

?>-i k-p|<2‘ mQ 


/T 3 


EE EL A q (u p ■ Vw p ')Wq dx 

q>— 1 p>g— 2 |p—p'[<l ^ 

= :/ 2 l + / 22 + / 2 3. 


Then, we estimate them as follows: 

|d 2 i| < V] V / |A 9 (u<p _ 2 • Vn>p)ui g | dx 
Jt 3 


P>Q k-p|<2 


<E E L |A 9 (u<q • Vwp)in g | dx - EE EL A q (u p > ■ Vw p )w q dx 

P>Q k-p|<2 ^ T3 p'>Q p>p'+2 |g-p|<2 ^ T3 

^II u <qIIooE E a pIHII 2 IKII 2 + X! Ilvlloo E E a pIKII 2 |KII 2 

P>Q \q— p|<2 p’>Q p>p' +2 |g—p|<2 


^ I|u<qIIoo X! A< ? Il^k1 !2 + ^2 A p'IK'lk E x iW w ' 

q>Q p’>Q q>p' 


2 

9 ll 2 • 
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Using definition of A Uyr and Jensen’s inequality, we obtain 

i-^2ii 1$ c r v y ' \q\q iinig 1I2 + c r v y ' x p i y ' 11 112 

q >Q p’>Q q>p' 

<c r v^X 2 q \\w q \\l + c r vY^ X 2 q \\w q \\l Y X q-p ' 
q>Q q>Q Q<p'<q 

< Cr-Z^ || VlD || 2 ■ 

We proceed estimating the second term. 

I-M < E E / |A 9 (up • Vw< p - 2 )w q \dx 

q>Q\q-p\<2 Jj3 

^ E E IKINK|| 2 Y VIK'll^ 

q>Q \q-p\<2 P'<q~2 


< C.~V Y^ X, 
q>Q 

< r ~ u Y x 

q>Q 


C r V l^Aq ' \\w q \\ 2 Y X l2 r Ill'll2 

p'<q-2 

C r V 2^W\ W q\\2 Y X p'\\ W P'h X Yp' 


p'<q -2 


< C r ^||VMj|| 2 - 

The last term J 2 3 is similar to 1 13 : 


|T 2 31 < Y, E / \ /X q( u P ® yjp)Vw q \dx 

q>Qp>q-2^ J3 

< E E L |A g (zi p <g> uip^Wql dx + Y, / |A g (wg ® wq)S7w q \ dx 

JT 3 JT 3 


p>Q Q<q<p+2 1 Q<q<Q-\-2' 

^ Y iKiWHih Y A <?iKii/^ + a q E H u qIIoo||wq||2||w,ii 2 

P>Q Q<q<P+ 2 Q<q<Q+2 


< 


C r vY X P "IHIk Y X q r \\w q \\2 + C r vX 2 Q E I|wq||2|K|| 2 


p>Q 


Q<q<p+2 


Q<q<Q +2 


<c r vY X p\\™ph Y x q\\ w qh x q-p + c r vX 2 Q E 


\w. 


q 112 


p>Q Q<«?<P+2 Q<q<Q+2 

< C r v\\X7w\\l. 

Therefore, we have 

(3-10) \I 2 \<c r v\\Vw\\ 2 . 

Combining (13.91) and (13.101) . we conclude that there exists an adimensional con¬ 
stant C > 0 that depends only on r, such that 

h +h < CctI/||Vm;||1 . 

Choosing c r := 1/(217) we infer from (13.81) that 


IK*)Hi <IH*o)lll - v [ \\Vw\\ 2 2 dT 

J to 


IMldr, 


'to 


<\\w(t 0 )\\l-u K 2 
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where ko = X' Thus 

IH*)ll! < lk(io)||ie—S(*-*°), 

which concludes the proof. 

□ 


4. Explicit estimates of the average determining wavenumber in 
terms of Kolmogorov’s dissipation wavenumber 


In this section we show that the average determining wavenumber has a uniform 
upper bound. Recall that 

A u ,r(t) := min{A q : X p r ||u P ||i,r < c r v , Mp> q and X~ ||u< g ||L=° < c r p, q £ N}, 

where r £ (2,3) and c r is an adimensional constant that depends only on r. To 
simplify the notations, denote A := A U)T and let Q be such that A q = A. First, it 
is easy to see that A(t) is locally integrablc. Indeed, whenever A(t) > Ao we have 
that one of the conditions in definition of A u ^ r is not satisfied, i.e., 

(4.11) \\uQ\\ r > c r vA x ~i, or ||m<q-i||oo > c^Xq-i = \c r vA. 

In the hrst case we use Bernstein’s inequality to obtain 

a^m^a^WuqWI 

< {c r v)- 2 A 2 \\u Q \\ 2 2 

< M- 2 \\Vu\\l 


In the second case, 


T<4(c r u)- 2 T- 1 ||u< Q _ 1 || 2 0 

< (c r i/) _2 yl _1 | Y IKII 

\q<Q 

< {CrV^A - 1 | ^ A|||u 9 || 2 

\q<Q 


= (crP) 2 Y x q\\ u qh X Q-q 
\q<Q 

< (c r v )~ 2 Y x >q\\l x Q l q 
q<Q 

<(c r u)- 2 ||Va|| 2 . 


We can now compare the average determining wavenumber with Kolmogorov’s dis¬ 
sipation wavenumber, often defined as 

* j:= (u) 4 ’ 

where e is the average energy dissipation rate 


(4.12) 


e-.= ^(\\Vu\\l) = ^- 3 


|| Vii(t)||| rfi, 
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and (•} is the time average. Then the above estimates on A yield the following 
bound: 


(4.13) 


(yl) — A 0 


< 


<l|Vu|||) 


= -L* = <jA, 


4\-3 


However, in this argument we did not take into account that the classical definition 
of Kolmogorov’s dissipation wavenumber assumes the absence of intermittency. In 
this case eddies occupy the whole region for each scale in the inertial range, and 
Bernstein’s inequalities (used to estimate A) are not sharp. On the other hand, 
the presence of intermittency requires the following modification of the classical 
definition (14.121) : 

(4.14) h vm H2>> 

where d £ [0,3] is the intermittency dimension. This parameter is chosen so that 


(4.15) 


<E A * 

q<Q 


-l+§-M(l-f), 


“qWr 


}<Ao 


’<E a : 

q<Q 


q 


)■ 


This can be done since by Bernstein’s inequality 

-1+7 11 


3-s -l+s,, 

Aq A q \\Uq\\2 


2 < Aq \\u q \\r < A“ || U, 


IWUqWl 


The intermittency dimension d measures a level of saturation of Bernsten’s inequal¬ 
ity (see Eli for similar definitions). The case d = 3 corresponds to Kolmogorov’s 
regime where at each scale in the inertial range eddies occupy the whole region. 
Note that —1 + d( 1 — -) < 0 since d < 3 and r < 3. Therefore, 


i i<Q 


(4.16) 


< £ A | 



Now, thanks to (14.111) . we have 

A -U-$+d(i- a)|| Ug|| 2 > {CrV) 2 A i+dii-l) 

or 

AA-^-^Wu^Wl > (c r v) 2 A^-t\ 
provided A > 1. Therefore, 

(c r i/) 2 yl 1+d ^ l¥ f ^ < H _1+ ^ +d ( 1_ ^||uQ|| 2 + 4H _1+d ( 1_ ^||u<Q_ 1 ||^ 0 , 
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whenever A > 1. Combining this with (14.1611 we obtain 
{CrvfA 1 ^-^ < A- 1+ -r +d ^\\u Q \\l + E A, 

q<Q 

= Y, A,- I+?+ ' <1 -'>ii«,iit 

q<Q 

Taking the time average and using (14.151) we arrive at 


- 1 + ' +ill -' ) ll<.X 




-l+fi+d(l-^) n ,, 2 


\q<Q 


>Aq 


<Ar f) E^iwii 

\q<Q / >Ao 


~ Aq 7 (HVwlll) 

2d f 

=: A ° ' ^ 

where (g) > \ 0 = (1 3 >a 0 ff)- Finally, using Jensen’s inequality, 


( y l) - Ao = (^}>A 0 = 


(4.17) 


< 




/ \ O / /i \ l - 1 - 1 — 77) \ 1+^(1— -p) 

{c r v Y(A) >Xo y ' 

(c r i/) 2 

^ M 2 J 

( £ \ 1 + <I(1-J) r+d(t-2) 

UsJ A o 


Comparing it with Kolmogorov’s dissipation wavenumber (14.141) we get 


w-ao s K r (i " ) Ao 7M . 

In case of extreme intermittency d = 0 the powers do not depend on r, so we can 
just choose r = 5/2 inferring that the average determining wavenumber is bounded 
by Kolmogorov’s dissipation wavenumber Kd '■ 

£ 

(yl) — Ao < —^ = Kd; when d = 0. 

However, in Kolmogorov’s regime d = 3, the bound becomes k^ + since r can only 
be taken less than 3. More precisely, for any d < 3 we have the following bound: 


(yl) — A 0 S C" 3 Aq 


— K d A o > 


where <5 = 2 r -3 can be arbitrarily small when r is chosen close to 3. 

5. Explicit estimates of the average determining wavenumber in 

TERMS OF THE GRASHOF NUMBER 

It is well known that the 3D Navier-Stokes equation possesses an absorbing ball 
in L 2 

B:={uGL 2 {T 3 ):\\u\\ 2 <R}, 
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where R is any number larger than 


and kq = 27tAq = 27 t/L, which can be 


expressed in terms of the adimensional Grashof number 


G := 


I H- 


2 V 2 


_i/o 

as R > i/k 0 ' G. Then for any Leray solution u(t) there exists to, depending only 
on 11 ( 0 ) 112 7 such that 

u(t) £ B Vt > t 0 . 

Then the evolutionary system consisting of Leray-Hopf solutions in the absorbing 
ball posseses a weak global attractor A, which has the following structure mm- 

A = {u(0) : u(-) is a complete bounded Leray-Hopf solution to the NSE}. 

The set A C B is the minimal L 2 weakly closed weakly attracting set, it is L 2 -weak 
omega limit of B (see [HO), and ||m ||2 < vk 0 A ' 2 G for all u £ A. 

Let u(t) be a trajectory on the global attractor A. To bound the average deter¬ 
mining wavenumber in terms of the Grashof number we use the energy inequality: 

rto~\~T pto~\~T 

0 < ||u(to + t)||| < limsup ||ix(t )111 —2v 

t—>to~\- J to 

ptQ-\-T 


< VK 0 l G 2 — V 


|| Vit(i)III dt 


|| Viz(t)||| dt + 2 

^ pto+T 


(/, u) dt 


to 


\ 2 h ~i dt. 


I to 


Therefore 


I 


rto+T /^2 

l|Vu(t)|||dt< — 


'to 


\ 2 H-1 


to 


0 


< +K0 V 2 G 2 . 

I KO 


Then gra imply 
(5.18) 


(A) - A 0 


< 


(llVulll) ^ G 


A ^ 


v* Tv 2 n o 


+ KqG 2 . 


To take into account intermittency, we can use (14.171) instead of (14.13[) to obtain 

D)-AoS(^)' 


l+d(l-£) 


2d 

\ r+d(r-2) 

A o 


'Aq (1 ^(llVnlll)' 


i+d(i-|> 


< KO 


G 2 


Tv 2 t 


+ G‘ 




where d £ [0,3] is the intermittency parameter from Section [I] In the case of 
extreme intermittency d = 0 this bound is the same as (15.181) . proportional to G 2 . 
However, in Kolmogorov’s regime where d = 3, the average determining number is 
bounded by G l+ . More precisely, 

4+5 


where <5 = ^ 


(H) — Ao <s ko 
+ 0 as r —> 3—. 


G 2 


Tv 2 h 


G“ 


d = 3, 
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6. The case of energy equality 

In this section we consider the case where one of the solutions of the 3D Navier- 
Stokes equations satisfies the energy equality, e.g., a steady state or a solution 
belonging to Onsager’s space L 3 ( 0, oo; B -\). 

Theorem 6.1. Let v(t) be a weak solution of the 3D NSE satisfying the energy 
equality, and Q(t) be such that A VtT (t) = X qo) for some r € (2,3). If u(t) is a 
Leray-Hopf weak solution such that 

(6-19) u(t)< Q{t) = v(t)< Q{t) , \/t > 0, 


then 


lim ||u(*) -v(t )\\ L 2 = 0. 


Proof. We know that v satisfies energy equality 

(6.20) f \-v 2 dx + v f f \Vv\ 2 dxdr = [ ^Vq dx ■ 
J t 3 2 J o Jj3 Jj3 2 


fv dx dr ; 


o Jt 3 


and u satisfies energy inequality 


(6.21) f \-u 2 dx + v f f \Vu\ 2 dxdr< f ^-u^dx+ f f fudxdr\ 
J T 3 2 J o J T 3 J T 3 ^ J o J T 3 

Computing the energy of the difference w := u — v 

f \-\w\ 2 dx + v f ( \X7w\ 2 dx dr 
J t 3 2 J o J T 3 

■j-u 2 dx + f \-v 2 dx + v f f |Vu| 2 dxdr + v l f |Vv| 2 da;dT 
2 J T 3 2 J q Jj 3 J Q Jf 3 


Jt 3 

uv dx — 2v 
It 3 


VuVi) dx dr, 


/o Jt 3 

combining (16.201) and (16.211) we conclude 

I ^-\w\ 2 dx + u f ( \X7w\ 2 dx dr 
J t 3 2 J o J T 3 

(6.22) < f i(ug + Uo)dx+ f ( f(u + v)dxdr 

J t 3 2 J o Jt 3 


/ 


uv dx — 2v 


VwV v dx dr. 


/T 3 JO J T 3 

As for the L 2 inner product of u and v, we have (see, e.g., [22] ) 


uqVq dx — uvdx — 2v 


(6.23) 


/T 3 


/ T 3 


VuVu dx dr 


o J T 3 


n w-Vvwdxdr— / / f(u + v)dxdr. 

' 3 Jo J T 3 


Combining (16.221) and (16.231) yields 


/T 3 


' dx + v f I | Vu>| 2 dx dr < f -|w(0)| 2 dx + f f \u> ■ Vvw\dx dr. 
Jo J T 3 J T 3 2 J o J T 3 
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It then follows from estimate (13.91) that for a small constant c r 


dx + v ( f |Viz;1 2 dxdr< f |w(0)| 2 da;. 

Jo J t 3 J T 3 

quality, we hav 

f \w\ 2 dx < f |w(0)| 

J T 3 J T 3 


J T 3 Jo JT 3 

Applying Poincare’s inequality, we have 


dx — kqv I / \w\ dxdr. 

0 J T 3 


Thus 


IKi)ll2<IK0)||^e 


2 —KQl/t 


□ 
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